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Abstract. We demonstrate the results of an attempt to match the two-soliton analytical 
solution with the numerically produced solutions of the Einstein field equations, that describe 
the spacetime exterior of rotating neutron stars, for arbitrary rotation. The matching procedure 
is performed by equating the first four multipole moments of the analytical solution to the 
multipole moments of the numerical one. We then argue that in order to check the effectiveness 
of the matching of the analytical with the numerical solution we should compare the metric 
components, the radius of the innermost stable circular orbit (Risco), the rotation frequency 
f2 = 4j| and the epicyclic frequencies Sl p , £1%. Finally we present some results of the comparison. 



1. Introduction 

Obtaining the geometry of the spacetime that surrounds a rotating distribution of mass in an 
exact form, is a problem of astrophysical relevance that remains unsolved. When modelling 
the astrophysical phenomena that we observe around rotating neutron stars, there are only a 
few choices we can make for the geometry we will use. We can either use one of the very 
well understood geometries of Schwarzschild or Kerr, which can only be used as very rough 
approximations, or we have to use a numerically obtained spacetime. 

Although there are various groups (see p], and for an extended list of numerical schemes see 
[2]), that have expertise in building relativistic models of astrophysical objects with adjustable 
physical characteristics and can construct the metric inside and outside such objects by solving 
numerically the full Einstein equations in stationary cases, having the metric of the spacetime 
in tabulated form with numerical values that correspond to the metric components at a selected 
grid is not always appropriate for studying the phenomena around the object. 

Thus, an analytical solution that could describe the geometry around such an object would 
be useful. Fortunately, nowadays, there is a large variety of analytical solutions of vacuum 
Einstein equations, which could be used as candidate metrics to describe well the exterior space- 
time of axisymmetric astrophysical objects. Ernst [3] formulated the Einstein equations in the 
case of axisymmetric stationary space-times long time ago, while Manko et al. and Sibgatullin 
IH 13 El El [H [9] have used various analytical methods to produce such spacetimes parameterized 
by various parameters that have a different physical context depending on the type of each 
solution. 

Such an analytical solution, appropriately matched to a neutron star, could be used to 
describe the stationary properties of the spacetime. That is, study the geodesies in the exterior 



of the neutron star. In other words, from the analytical solution, we could obtain bounds of 
motion for test particles, the location of the innermost stable circular orbit (ISCO), the rotation 
frequency of the circular orbits on the equatorial plane as well as the epicyclic frequencies around 
them. 

These properties that are of astrophysical relevance and characteristic of a spacetime, will 
be used to compare the proposed analytical solution to the numerically produced solutions that 
describe the spacetime exterior to a neutron star. 

The rest of the paper is organized as follows: In Sec. [5] the proposed analytical solution (two- 
soliton) is briefly presented and some of its properties analyzed. In Sec. Owe present at first the 
matching criteria and then the comparison criteria and some results. 



2. An Analytical axially symmetric solution: The Two Soliton solution. 

The vacuum region of a stationary and axially symmetric spacetime can be described by the 
Papapetrou line element [10] 



ds 2 = -f {dt - ujd^) 2 + r l [e 27 (dp 2 + dz 2 ) + p 2 d(/) 2 ] , (1) 

were the functions /, oj and 7 are functions of the Weyl-Papapetrou coordinates (p, z). Using 
this metric, the Einstein field equations in vacuum reduce to the Ernst equation |3| 

Re(£)V 2 £ = V£ • V£ , (2) 

were the Ernst potential £, is a complex function of the metric functions. 

A general procedure for generating solutions of the Ernst equation was developed by 
Sibgatullin and Manko El El Hj. The solution of the Ernst equation is produced from a 
choice of the Ernst potential on the axis of the form 

£(p = 0,z) = e(z) = ^, (3) 

were the functions P(z), R(z) are polynomials of order n with complex coefficients in general. 
In [UJ, one can find a detailed description of the procedure for generating the solutions using 
the parameters of the Ernst potential. It is shown in [BJ that the metric functions are finally 
expressed in terms of some determinants, which depend on the parameters £ n , which are the 2n 
roots of the equation e(z) + e(z)* = 0, the parameters e n , f n and f3 n , which are defined from 
the expressions 



n n p 



k=l MK fc=l 



and finally the expressions R n = \J p 2 + (z — £ n ) 2 , which are functions of the coordinates (p, z). 

The vacuum two-soliton solution (proposed by Manko 0]) is a special case of the previous 
general axisymmetric solution that is obtained from the ansatz (see also |11| ) 



(z — M — ia)(z + ib) — k 
6{ - Z > = (z + M -ia)(z + ib)-k ( ' 

where all the parameters are real while the parameters M, a are the mass and the reduced 
angular momentum -j^ respectively. The first five mass and angular momentum moments of the 
corresponding spacetime are: 



M = M, Mi = 0, M 2 = -(a 2 - k)M, M 3 = 0, 
M 4 = (a A -(3a 2 -2ab + b 2 )k + k 2 + ^(2a 2 -k)M 2 ^jM 
J = 0, Ji = aM, J 2 = 0, 

J 3 = -(a 3 - (2a-b)k)M, J 4 = 0. (5) 
For the ansatz ([!]), the equation e(z) + e(z)* = corresponds to: 

z 4 - (M 2 -a 2 -b 2 + 2k)z 2 + (k- abf - b 2 M 2 = 0, (6) 

where the coefficients take real values. That means that the roots can either be real, or conjugate 
pairs. We can also see that the four roots of ([6]) are of the form 

£i = -£3 = £+, 6 = = 

Defining the parameters d = y/(k — ab) 2 — b 2 M' 2 and k± = V M 2 — a 2 — b 2 + 2k ± 2d, we can 
express the roots of (0) as 

£± = ^(k+±«_). (7) 
Using these parameters we redefine R n as 

R ± = ^p 2 + {z±t + ) 2 , r± = yW(^±£-) 2 - (8) 

Depending on the choice of the parameters used in (jj]), there can be obtained various types 
of solutions. For example, if one sets the value of k equal to zero, for arbitrary values of the 
parameter 6, the two-soliton corresponds to the Kerr metric of mass M and reduced angular 
momentum a (if \a\ < M) or the hyperextreme Kerr (if \a\ > M). Further more, if one chooses 
a equal to zero, then the two-soliton corresponds to the Schwarzschild solution of mass M . 

There is also another interesting special case for the two-soliton solution. If the parameter d 
is real and we have the conditions M 2 - a 2 - b 2 + 2k + 2d > and M 2 - a 2 - b 2 + 2k - 2d = 0, 
then the parameter n+ is real and the parameter k_ = 0. In that case there is a degeneracy and 
there are only two real roots £1^2 = ±£± = i^f - This solution corresponds to the Manko et al. 
solution [5l [7] . 

In brief, the two-soliton solution can produce a very rich family of analytical solutions 
including the classical solutions of Schwarzschild and Kerr as well as the Manko et al. solution 
O [7j that has been used in [12[ [T3] to approximately match the exterior of a rotating neutron 
star. 

At this point we can point out the advantages of the two-soliton solution against the Manko 
et al. solution [5j [7]. If we assume a constant mass M, the two-soliton solution has a three 
dimensional parameter space (a, b, k). In that space the plane k = includes all Kerr solutions 
of mass M while the b axis corresponds to the Schwarzschild solution of mass M. 

As we have mentioned earlier, the Manko et al. solution [21 E] is obtained if we impose the 
constrain M 2 — a 2 — b 2 + 2k — 2d = to the parameters. This is analogous to having a two 
dimensional surface in the parameter space. There are two things that could be pointed out 
about that surface, which is drawn in Fig. [TJ The first is that the axis of b (the blue line) that 
corresponds to the Schwarzschild solution doesn't cross the surface and the second is that the 
Kerr solutions with M > \a\, that lie in the parameter space on the plane k = between the 
two red lines parallel to the b axis, do not intersect the surface either. 



Figure 1. The first two plots show the parameter space (a, b, k) of the two-soliton solution for 
a specific mass (M = 1) from two different view points. The 2D surface plotted is the constrain 
of k which corresponds to the Manko et al. solution [51 [7] used by Berti and Stergioulas [12] . 
The third plot is the k = plane of the parameter space which corresponds to all Kerr solutions 
of the two-soliton solution. The plots clearly show that the Manko et al. solution has no set of 
parameters that describe the case of k = which corresponds to the Kerr and the Schwarzschild 
(for a = 0) solutions, since there is no intersection between the constrain of k and the Kerr plane 
within the appropriate range of parameters. The two hyperbola plotted on the plane k = are 
the points were the constrain of k touches the plane in a tangent manner. As we can see these 
hyperbola correspond to |o| > M. 

Thus, this solution can reduce to Schwarzschild or Kerr only formally. These are important 
drawbacks and were pointed out by Berti and Stergioulas [12] . Particularly important is the 
fact that it can only be matched with rapidly rotating models of neutron stars, since in the limit 
of zero rotation it has a non- vanishing quadrupole moment Q = — ^ ^ngg~^ • These problems 
do not exist for the two-soliton solution, that can represent neutron stars of arbitrary rotation. 
Another important aspect of the two-soliton is that, as can be seen from the multipole moments 
([5]), the solution can describe spacetimes with an arbitrary quadrupole M2 and an arbitrary 
current octupole J3, since every parameter is introduced linearly in the first moment it appears 
and thus it can always be determined by the first four non-zero moments. 

3. Matching and comparison of the analytical to the numerical metric. 

3.1. Matching of the analytical to the numerical solution. 

In matching the analytical solution to the numerical one it is desirable to find a criterion that 
will be characteristic of the whole structure of the spacetime and not of a small region of it. 
That is, the matching conditions should be global and not local. Berti and Stergioulas [12] have 
argued that these global conditions should be the matching of the multipole moments. We will 
work under the same assumption. The multipole moments of a stationary axially symmetric 
spacetime are connected to the Ernst potential on the axis of symmetry and the spectrum of 
the moments fully specify it. On the other hand, when the Ernst potential on the axis is given, 
there exists a spacetime which is fully specified by that Ernst potential [14]. Thus, the multipole 
moments are characteristic of a spacetime and can be used as a global matching condition. 

When a numerical spacetime for a neutron star model is constructed from a numerical 
algorithm, it is possible to evaluate its mass moments M, Q, ... and angular momentum moments 
J, S3, ... (for further discussion see Berti and Stergioulas |12j). These numerically evaluated 
moments can be used to impose the matching conditions to the analytical spacetime. The first 



four nonzero multipole moments of the two-soliton solution are: 

M = M , Ji = aM, 
M 2 = -{a 2 -k)M , J 3 = -(a 3 - (2a-b)k)M, (9) 

where we can see that once we specify the mass and the angular momentum of the spacetime, 
the parameter k is uniquely determined by the quadrupole moment Q = M<i and the parameter 
b is uniquely determined by the current octupole S3 = J3. Thus, having specified the parameters 
of the two-soliton, we have completely determined the analytical spacetime that can be used to 
describe the exterior of the neutron star. What remains to be seen is how well do the properties 
of the analytical spacetime compare to those of the numerical one. 



3.2. Criteria for the comparison of the analytical to the numerical spacetime. 
To compare the analytical to the numerical spacetime, one should use criteria that are 
characteristic of the geometric structure of the spacetime. The criteria should also be related to 
properties of the spacetime that could potentially be measured in astrophysical phenomena. A 
good fit of such properties would mean that they could be used to identify a spacetime through 
observational data and then associate it with properties of the compact object that generates it. 

The first and more obvious criterion one can use is how well the respective components of 
the analytical and the numerical metric match. If the relative difference of the analytical metric 
components 

9tt = ~f, 9t<t> = M 3<H> = /~V ~ f^ 2 ( 10 ) 
to the corresponding numerical is small, one could consider the former as a good substitute of 
the latter. The comparison of the component can be performed through the circumferential 
radius R c irc = ^/9~H> which is as an indicator of how well the measurements of the circumferences 
in the corresponding spaces (analytical and numerical) with a specific value of the coordinate p, 
compare with each other. 

The second criterion is the innermost stable circular orbit (ISCO). Test particles that move 
on geodesies on the equator follow trajectories that are given by the equation 



9pp[^z) =1 "a- 2 = V(p) (11) 



'dpy _ E 2 g H + 2ELg t<p + L 2 g tt 

where E and L are the conserved energy per unit mass and angular momentum (parallel to 
the axis component) per unit mass and V(p) is the effective potential. Since we are interested 
in circular orbits, the conditions for such orbits are dp/dr = and d 2 p/dr 2 = which are 
equivalent from (jlip to the conditions for the effective potential, V{p) = 0, dV{p)/dp = 0. 
The ISCO is determined if we additionally impose the condition that the radius p of the circular 
orbit is also a turning point of the potential, that is d 2 V(p)/dp 2 = 0. The location of the ISCO 
is of astrophysical importance since it indicates the radius at which matter orbiting a compact 
object can no longer maintain its orbit and plunges into the object. That radius would be the 
internal radius of accretion disks formed around such a compact object. 

The third criterion is the rotation frequency of circular orbits on the equatorial plane. The 
rotation frequency Q = -£ is determined from (jlip and the conditions for circular orbits, 
V{p) = 0, dV{p)/dp = 0, and is given by the equation 



-9t<t>,p+ \/W,p) -9tt, P 9u,P /10 , 
il(p) = . (12) 

94><t>,p 



The final criterion are the epicyclic frequencies Q p , £l z . These are the frequencies of precession 
of the periastron and the orbital plane of a perturbed circular orbit, respectively. They can be 
evaluated by perturbing the equations of motion around a circular orbit 



where V(p, z) is the effective potential defined in (jlip maintaining its z dependence. The final 
expressions for the epicyclic frequencies are 



where a is p or z and the expression is evaluated at z = and around circular orbits of p = pq. 

These criteria can be related to astrophysical phenomena and in particular, they can be re- 
lated to properties and observables of accretion disks. One such observable is the quasi-periodic 
oscillations (QPOs) with modulations in the range of kHz that have been observed in the X-ray 
flux from accreting X-ray pulsars [15]. The association of QPOs with the above frequencies 
could be used to determine the properties of the central object from the values of the various 
frequencies. These properties such as mass, angular momentum and mass quadrupole, could 
also be used to gain insight on the structure of the central object. Thus, an analytical spacetime 
that fits well all these parameters could be practically considered equivalent to the numerical 
one and a useful tool for doing astrophysics. 



3. 3. Results of the comparison. 

In the comparison of the metric components, we have used the two corresponding Manko et al. 
solutions that were used by Berti and Stergioulas [12] , as a reference. The comparison of the 
metrics is performed on the equatorial plane and on the axis of rotation. Here we will present 
only one rotating neutron star model, since these results are typical for all models. The chosen 
example is the one in table 3 of the Berti and Stergioulas |12] work, that corresponds to the 
model with the fastest rotation of the sequence of maximum mass in the non-rotating limit, for 
the FPS equations of state (EOS). 

In Fig. [2] we have plotted the relative difference of the metric components gu and <?m and 
of the circumferential radius R c irc a s a function of coordinate radius p on the equatorial plane 
(the T> operator is defined as T> = (numerical - analytical) /numerical). In the first two graphs 
we have used for comparison the Two Soliton solution and the two corresponding solutions of 
Manko et al. which are distinguished by their value of the parameter b of the Ernst potential. 
We remind here that the Ernst potential of the Manko et al. solution is given by the expression 




(13) 




(14) 




(z - M - ia)(z + ib) + d - 5 - ab 



(15) 



(z + M - ia)(z + ib) + d - 5 - ab 





Figure 2. The first two plots show the rel- 
ative difference of the three analytical met- 
rics with the numerical for the components 
gu and g t ^ on the equatorial plane. The an- 
alytical metrics used are the two-soliton and 
the two Manko et al. [3E] that correspond to 
the two values of the parameter b (see Berti 
and Stergioulas [12 J). The third plot shows 
the relative difference of the circumferential 
radius on the equatorial plane for the two- 
soliton and the Manko et al. with the negative 
value of b (the two solutions almost coincide). 
As we can see, the errors are well under 1 per- 
cent (for the green line) on the surface and fall 
outwards. The oscillations that appear in the 
second plot are due to the behavior of the nu- 
merical solution. 
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where there are only three independent parameters M, a, b and for a particular quadrupole 
moment of the numerical solution (with fixed M, a) there may be no corresponding value of b 
or there may be two values. For these two values of b we can calculate two metrics which are the 
ones shown in the graphs. As the second graph shows, the one of the two Manko et al. solutions 
does not fit well with the numerical metric and thus it describes it poorly [12]. On the other 
hand, the other solution (corresponding to the negative value of b) is as good as the two-soliton 
metric. Both metrics are better than 1 percent on the equator and the relative difference falls 
with increasing radial coordinate p. The relative difference shows a kink at about 100km which 
is an artifact due to a crossing of the numerical and the analytical solutions and has no further 
physical meaning. In the third graph of Fig. [2] we have used only the one of the Manko et al. 
metrics (the one with negative b) and the two-soliton. The relative difference for the R c %rc is 
under 0.1 percent on the surface and falls with increasing radial coordinate p. 

In Fig. [3]we have plotted the relative difference of gu on the axis of symmetry. The plot shows 
the same good behavior of the analytical metric on the axis as on the equatorial plane. The 
behavior on the axis of the analytical metric was a matter of confusion that has been clarified 
in 




Figure 3. The plot shows the relative 
difference of the three analytical metrics with 
the numerical one for the component gu on 
the axis of symmetry. The relative difference 
is under 1 percent on the surface. The 
behavior is similar to the behavior on the 



logiz) (km) 

equatorial plane. 



As we can see, both the Manko et al. and the two-soliton fit well the numerical metric. But 
for the Manko et al. there are constrains. If the rotation rate of the neutron star becomes 
smaller than a critical value (for typical values of the rotation parameter j = J/M 2 see j!2|). 
there cannot be found any solutions (based on the first four multipole moments) of the type 
of Manko et al. that could match to the exterior of the neutron star. On the other hand, the 
two-soliton can describe arbitrary rotations as was pointed out and the good fit to the numerical 
metric extends to all rotations. 
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Figure 4. Relative difference of the 
analytical with the numerical ISCO for the 
sequence of maximum mass in the non- 
rotating limit for the equation of state FPS. 
The maximum difference of the co-rotating 
ISCO is 5.7 percent for the fastest rotating 
model that has a co-rotating ISCO and the 
maximum difference for the counter rotating 
ISCO is less than 1 percent for the fastest 
rotating model. 



For the comparison of the ISCO we have used a sequence of models of neutron stars with 
various rotations. The sequence is the one in table 3 of Berti and Stergioulas [12] that corresponds 
to the model of maximum mass in the non-rotating limit for the FPS equation of state. The 
comparison of the Risco is shown in Fig. [H where we have plotted the relative difference of 
the numerical ISCO with the analytical one as a function of the rotation parameter j = J/M 2 . 
The negative values of j correspond to the counter-rotating ISCOs which are located further 
from the star surface. In some cases the co-rotating ISCO is under the surface of the star so, 
there is no numerical value for it, while that never happens for the counter-rotating ISCO. The 
co-rotating analytical ISCO is systematically less accurate than the counter-rotating one and 
that is to be expected since the latter is located further from the star. The difference in the 
ISCOs is less than 5.7 percent for the co-rotating and less than 1 percent for the counter-rotating 



orbits. The same general results are true for the Manko et al. metric for rotation parameters 
above the critical value. We should point out that the results of Fig. [3] are different than the 
ones shown in Fig. 8 of [12] . The difference is on the counter-rotating branch and it is due to 
a calculation slip that resulted to overestimation of the error in the counter-rotating ISCOs in 
|12| . For example the fastest rotating model of fig. 7 in [12j is given an error of 10 percent while 
the true error is of the order of 1.2 percent. Thus the ISCO comparison is more favorable than 
it was originally estimated. 

The next property that we compare is the rotation frequency of the circular orbits. The 
relative difference of f2 on the equatorial plane is plotted as a function of the coordinate distance 
p and is shown in Fig. [SJ The relative difference is less than 1 percent on the surface of the 
star and reduces with increasing distance p. In this model, since the ISCO is located under the 
surface of the neutron star, there is no problem to evaluate f2 from the surface outwards. The 
analytically evaluated f2, as Fig. [5] shows, fits well the numerical and that conclusion also 
holds for all rotations. 




Figure 5. Relative difference in Vt of the two- 
soliton with the numerical metric. As in the 
plots of the various metric components, the 
relative difference is under 1 percent on the 



^gp(km) surface and falls outwards. 



3.4- Conclusions 

That was a brief presentation of our results on comparing the analytical two-soliton metric with 
numerical metrics describing the exterior of neutron stars. We plan to present a more thorough 
comparison between the two kinds of metrics in a forthcoming paper [T7] where the comparison 
of the epicyclic frequencies will be demonstrated as well. 

The results presented here suggest that the two-soliton metric is a very good candidate 
for describing the exterior spacetime of a neutron star of arbitrary rotation. Therefore, we 
could use this metric to gain information about the neutron star through observables related to 
astrophysical processes in the vicinity of the star. 
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